Relating target fish DNA concentration to community composition analysis in freshwater fish via metabarcoding.

Metabarcoding has been widely accepted as a useful tool for biodiversity assessment based on eDNA. The method allows for the detection of entire groups of organisms in a single sample, making it particularly applicable in aquatic habitats. The high sensitivity of the molecular approaches is especially beneficial in detecting elusive and rare fish species, improving biodiversity assessments. Numerous biotic and abiotic factors that affect the persistence and availability of fish DNA in surface waters and therefore affecting species detectability, have been identified. However, little is known about the relationship between the total fish DNA concentration and the detectability of differential abundant species. In this study three controlled mock-community DNA samples (56 individual samples) were analyzed by (i) metabarcoding (MiSeq) of 12S rDNA (175 bp) and by (ii) total freshwater fish DNA quantification (via qPCR of 12S rDNA). We show that the fish DNA quantity affects the relative abundance of species-specific sequences and the detectability of rare species. In particular we found that samples with a concentration between 1000 pg/µL down to 10 pg/µL of total fish DNA revealed a stable relative frequency of DNA sequences obtained for a specific fish species, as well as a low variability between replicates. Additionally, we observed that even in complex mock-community DNA samples, a total fish DNA concentration of 23 pg/µL was sufficient to reliably detect all species in every replicate, including three rare species with proportions of ≤0.5 %. We also found that the DNA barcode similarity between species can affect detectability, if evenness is low. Our data suggest that the total DNA concentration of fish is an important factor to consider when analyzing and interpreting relative sequence abundance data. Therefore, the workflow proposed here will contribute to an ecologically and economically efficient application of metabarcoding in fish biodiversity assessment.

Silver nanoparticles adversely affect the swimming behavior of European Whitefish (Coregonus lavaretus) larvae within the low µg/L range.

The aim of this study was to determine the effects of silver nanoparticles (AgNPs; speciation: NM-300 K) in the lab on the behavior of larvae in European Whitefish (Coregonus lavaretus), a relevant model species for temperate aquatic environments during alternating light and darkness phases. The behavioral parameters measured included activity, turning rate, and distance moved. C. lavaretus were exposed to AgNP at nominal concentrations of 0, 5, 15, 45, 135, or 405 µg/L (n = 33, each) and behavior was recorded using a custom-built tracking system equipped with light sources that reliably simulate light and darkness. The observed behavior was analyzed using generalized linear mixed models, which enabled reliable detection of AgNP-related movement patterns at 10-fold higher sensitivity compared to recently reported standard toxicological studies. Exposure to 45 µg/L AgNPs significantly resulted in hyperactive response patterns for both activity and turning rates after an illumination change from light to darkness suggesting that exposure to this compound triggered escape mechanisms and disorientation-like behaviors in C. lavaretus fish larvae. Even at 5 µg/L AgNPs some behavioral effects were detected, but further tests are required to assess their ecological relevance. Further, the behavior of fish larvae exposed to 135 µg/L AgNPs was comparable to the control for all test parameters, suggesting a triphasic dose response pattern. Data demonstrated the potential of combining generalized linear mixed models with behavioral investigations to detect adverse effects on aquatic species that might be overlooked using standard toxicological tests.

A comprehensive DNA barcode inventory of Austria's fish species.

Austria is inhabited by more than 80 species of native and non-native freshwater fishes. Despite considerable knowledge about Austrian fish species, the latest Red List of threatened species dates back 15 years and a systematic genetic inventory of Austria's fish species does not exist. To fulfill this deficit, we employed DNA barcoding to generate an up-to-date and comprehensive genetic reference database for Austrian fish species. In total, 639 newly generated cytochrome c oxidase subunit 1 (COI) sequences were added to the 377 existing records from the BOLD data base, to compile a near complete reference dataset. Standard sequence similarity analyses resulted in 83 distinct clusters almost perfectly reflecting the expected number of species in Austria. Mean intraspecific distances of 0.22% were significantly lower than distances to closest relatives, resulting in a pronounced barcoding gap and unique Barcode Index Numbers (BINs) for most of the species. Four cases of BIN sharing were detected, pointing to hybridization and/or recent divergence, whereas in Phoxinus spp., Gobio spp. and Barbatula barbatula intraspecific splits, multiple BINs and consequently cryptic diversity were observed. The overall high identification success and clear genetic separation of most of the species confirms the applicability and accuracy of genetic methods for bio-surveillance. Furthermore, the new DNA barcoding data pinpoints cases of taxonomic uncertainty, which need to be addressed in further detail, to more precisely assort genetic lineages and their local distribution ranges in a new National Red-List.

Incubation media modify silver nanoparticle toxicity for whitefish (Coregonus lavaretus) and roach (Rutilus rutilus) embryos.

Toxicological studies were performed to examine silver nanoparticle (AgNP, size: 14.4 ± 2.5 nm) transformation within three different test media and consequent effects on embryos of whitefish (Coregonus lavaretus) and roach (Rutilus rutilus). The test media, namely ASTM very hard water, ISO standard dilution medium, and natural lake water differed predominantly in ionic strength. Total silver was determined using inductively coupled plasma mass spectrometry (ICP-MS), while AgNPs were characterized by transmission electron microscopy and single particle ICP-MS. Silver species distributions were estimated via thermodynamic speciation calculations. Data demonstrated that increased AgNP dissolution accompanied by decreasing ionic strength of the test medium did not occur as noted in other studies. Further, other physicochemical parameters including AgNP size and metallic species distribution did not markedly affect AgNP-induced toxicity. Irrespective of the test medium, C. lavaretus were more sensitive to AgNP exposure (median lethal concentration after 8 weeks: 0.51-0.73 mg/L) compared to R. rutilus, where adverse effects were only observed at 5 mg/L in natural lake water. In addition, AgNP-induced toxicity was lower in the two standard test media compared to natural lake water. Currently, there are no apparent studies assessing simultaneously the sensitivity of C. lavaretus and R. rutilus to AgNP exposure. Therefore, the aim of this study was to (1) investigate AgNP-induced toxicity in C. lavaretus and R. rutilus cohabiting in the same aquatic environment and (2) the role played by test media in the observed effects of AgNPs on these aquatic species.

DNA Transposon Expansion is Associated with Genome Size Increase in Mudminnows.

Genome sizes of eukaryotic organisms vary substantially, with whole-genome duplications (WGD) and transposable element expansion acting as main drivers for rapid genome size increase. The two North American mudminnows, Umbra limi and Umbra pygmaea, feature genomes about twice the size of their sister lineage Esocidae (e.g., pikes and pickerels). However, it is unknown whether all Umbra species share this genome expansion and which causal mechanisms drive this expansion. Using flow cytometry, we find that the genome of the European mudminnow is expanded similarly to both North American species, ranging between 4.5 and 5.4 pg per diploid nucleus. Observed blocks of interstitially located telomeric repeats in U. limi suggest frequent Robertsonian rearrangements in its history. Comparative analyses of transcriptome and genome assemblies show that the genome expansion in Umbra is driven by the expansion of DNA transposon and unclassified repeat sequences without WGD. Furthermore, we find a substantial ongoing expansion of repeat sequences in the Alaska blackfish Dallia pectoralis, the closest relative to the family Umbridae, which might mark the beginning of a similar genome expansion. Our study suggests that the genome expansion in mudminnows, driven mainly by transposon expansion, but not WGD, occurred before the separation into the American and European lineage.

Complex and divergent histories gave rise to genome-wide divergence patterns amongst European whitefish (Coregonus lavaretus).

Pleistocene glaciations dramatically affected species distribution in regions that were impacted by ice cover and subsequent postglacial range expansion impacted contemporary biodiversity in complex ways. The European whitefish, Coregonus lavaretus, is a widely distributed salmonid fish species on mainland Europe, but in Britain it has only seven native populations, all of which are found on the western extremes of the island. The origins and colonization routes of the species into Britain are unknown but likely contributed to contemporary genetic patterns and regional uniqueness. Here, we used up to 25,751 genome-wide polymorphic loci to reconstruct the history and to discern the demographic and evolutionary forces underpinning divergence between British populations. Overall, we found lower genetic diversity in Scottish populations but high differentiation (FST  = 0.433-0.712) from the English/Welsh and other European populations. Differentiation was elevated genome-wide rather than in particular genomic regions. Demographic modelling supported a postglacial colonization into western Scotland from northern refugia and a separate colonization route for the English/Welsh populations from southern refugia, with these two groups having been separated for more than ca. 50 Ky. We found cyto-nuclear discordance at a European scale, with the Scottish populations clustering closely with Baltic population in the mtDNA analysis but not in the nuclear data, and with the Norwegian and Alpine populations displaying the same mtDNA haplotype but being distantly related in the nuclear tree. These findings suggest that neutral processes, primarily drift and regionally distinct pre-glacial evolutionary histories, are important drivers of genomic divergence in British populations of European whitefish. This sheds new light on the establishment of the native British freshwater fauna after the last glacial maximum.

Monitoring spawning migrations of potamodromous fish species via eDNA.

Potamodromous fish are considered important indicators of habitat connectivity in freshwater ecosystems, but they are globally threatened by anthropogenic impacts. Hence, non-invasive techniques are necessary for monitoring during spawning migrations. The use of environmental DNA (eDNA) potentially facilitates these efforts, albeit quantitative examinations of spawning migrations remain so far mostly uncharted. Here, we investigated spawning migrations of Danube bleak, Alburnus mento, and Vimba bream, Vimba vimba, and found a strong correlation between daily visual fish counts and downstream eDNA signals obtained from filtered water samples analysed with digital PCR and end-point PCR coupled with capillary electrophoresis. By accounting for daily discharge fluctuations, it was possible to predict eDNA signal strength from the number of migrating fish: first, the whole spawning reach was taken into account. Second, the model was validated using eDNA signals and fish counts obtained from the upper half of the examined river stretch. Consequently, fish counts and their day-to-day changes could be described via an eDNA-based time series model for the whole migration period. Our findings highlight the capability of eDNA beyond delivering simple presence/absence data towards efficient and informative monitoring of highly dynamic aquatic processes such as spawning migrations of potamodromous fish species.

Zebrafish larvae show negative phototaxis to near-infrared light.

Zebrafish larvae (Danio rerio) are among the most used model species to test biological effects of different substances in biomedical research, neuroscience and ecotoxicology. Most tests are based on changes in swimming activity of zebrafish larvae by using commercially available high-throughput screening systems. These systems record and analyse behaviour patterns using visible (VIS) and near-infrared (NIR) light sources, to simulate day (VIS) and night (NIR) phases, which allow continuous recording of the behaviour using a NIR sensitive camera. So far, however, the sensitivity of zebrafish larvae to NIR has never been tested experimentally, although being a critical piece of information for interpreting their behaviour under experimental conditions. Here, we investigated the swimming activity of 96 hpf (hours post fertilization) and 120 hpf zebrafish larvae under light sources of NIR at 860 nm and at 960 nm wavelength and under VIS light. A thermal source was simultaneously presented opposite to one of the light sources as control. We found that zebrafish larvae of both larval stages showed a clear negative phototactic response towards 860 nm NIR light and to VIS light, but not to 960 nm NIR light. Our results demonstrated that zebrafish larvae are able to perceive NIR at 860 nm, which is almost identical to the most commonly used light source in commercial screening systems (NIR at 850 nm) to create a dark environment. These tests, however, are not performed in the dark from the zebrafish´s point of view. We recommend testing sensitivity of the used test organism before assuming no interaction with the applied light source of commonly used biosensor test systems. Previous studies on biological effects of substances to zebrafish larvae should be interpreted with caution.

Thermal experience during embryogenesis contributes to the induction of dwarfism in whitefish Coregonus lavaretus.

Ecotype pairs provide well-suited model systems for study of intraspecific phenotypical diversification of animals. However, little is still known about the processes that account for the development of different forms and sizes within a species, particularly in teleosts. Here, embryos of a normal-growing 'large' form and a dwarf form of whitefish Coregonus lavaretus were incubated at two temperatures that are usually experienced at their own spawning sites (2°C for the normal and 6°C for the dwarf form). All fish were subjected to similar thermal treatment after hatching. The present data demonstrate for the first time that different thermal experience in embryonic life has lasting effects on body and muscle growth of this ecotype pair and contributes to the development of the dwarf form. Thus, juvenile fish of the regular form are much smaller and have less muscle mass when pre-hatching thermal conditions were similar to those typical for the spawning sites of the dwarf form (6°C) than when subjected to conditions of their own spawning sites (2°C). Surprisingly, fish of the dwarf form exhibit a similar pattern of response to thermal history (2°-fish much larger than 6°-fish), indicating that in their case, normal spawning site temperature (6°C) is indeed likely to act as a growth limiting factor. Results also demonstrate that the hypertrophic and hyperplastic muscle growth modes are similarly affected by thermal history. Immunolabelling experiments for Pax7, H3P and Mef2 provide evidence that the cellular mechanisms behind the increased growth rates after cold incubation in both ecotypes are increased proliferation and reduced differentiation rates of muscle precursor cells. This is of major significance to aspects of ecological and developmental biology and from the evolutionary perspective.

Morphological divergence of lake and stream Phoxinus of Northern Italy and the Danube basin based on geometric morphometric analysis.

Minnows of the genus Phoxinus are promising candidates to investigate adaptive divergence, as they inhabit both still and running waters of a variety of altitudes and climatic zones in Europe. We used landmark-based geometric morphometric methods to quantify the level of morphological variability in Phoxinus populations from streams and lakes of Northern Italy and the Danube basin. We analyzed body shape differences of populations in the dorsal, lateral, and ventral planes, using a large array of landmarks and semilandmarks. As the species identification of Phoxinus on morphological characters is ambiguous, we used two mitochondrial genes to determine the genetic background of the samples and to ensure we are comparing homogenous groups. We have found significant body shape differences between habitats: Minnow populations inhabiting streams had a deeper body and caudal peduncle and more laterally inserted pectoral fins than minnows inhabiting lakes. We have also found significant body shape differences between genetic groups: Italian minnows had deeper bodies, deeper and shorter caudal peduncles, and a shorter and wider gape than both groups from the Danube. Our results show that the morphology of Phoxinus is highly influenced by habitat and that body shape variation between habitats was within the same range as between genetic groups. These morphological differences are possibly linked to different modes of swimming and foraging in the respective habitats and are likely results of phenotypic plasticity. However, differences in shape and interlandmark distances between the groups suggest that some (though few) morphometric characters might be useful for separating Phoxinus species.

Hybridization and restricted gene flow between native and introduced stocks of Alpine whitefish (Coregonus sp.) across multiple environments.

Translocations of Baltic whitefish (Coregonus sp.) into Austrian Alpine lakes have created 'artificial hybrid zones', threatening the genetic integrity of native lineages. We evaluate the genetic structure of Coregonus in Austrian lakes and characterize hybridization and introgression between native and introduced lineages. Fifteen populations (N=747) were assessed for allelic variation at eight microsatellite loci and a reduced set (N=253) for variation across two mtDNA genes (cyt b and NADH-3). Bayesian approaches were used to estimate individual admixture proportions (q-values) and classify genotypes as native, introduced or hybrids. q-value distributions varied among populations highlighting differential hybridization and introgression histories. Many lakes revealed a clear distinction between native and introduced genotypes despite hybridization, whereas some locations revealed hybrid swarms. Genetic structure among lakes was congruent with morphological divergence and novelty raising speculation of multiple taxa, including a population south of the Alps, outside the putative native range of Coregonus. Although statistically congruent with inferences based on nuclear markers, mitochondrial haplotype data was not diagnostic with respect to native and non-native lineages, supporting that the Alpine region was colonized post-glacially by an admixture of mtDNA lineages, which coalesce >1 Ma. Mechanisms promoting or eroding lineage isolation are discussed, as well as a high potential to conserve native Alpine lineages despite the extensive historical use of introduced Baltic stocks.

Changing handling times during feeding and consequences for prey size selection of 0+ zooplanktivorous fish.

The interrelationship of fish size, prey size and handling time within a 15-min feeding period was studied in three size groups of 0 + roach, Rutilus rutilus, and bleak, Alburnus alburnus. Four size classes of cladoceran prey were used to measure changes in feeding rate and handling time from initial rapid feeding to sustained feeding. Observed differences in increase of handling time between prey size classes led to a change in the prey profitability ranking of those size classes within the first 2 min of the experiments. A 1-min feeding period is interpreted as reflecting an intermediate motivational status between extreme hunger and satiation. The use of average handling times for this period revealed a substantial change in prey profitability estimates compared to previous studies which used handling times based on short-term (a few seconds up to 1 min) feeding. It is not the largest prey items a fish can handle and swallow that are most profitable, but prey of intermediate size. By this approach a closer fit between expectations derived from optimal foraging theory and empirical data on prey size selection of 0 + zooplanktivorous fish is qualitatively achieved. Optimal prey size was found to be close the mouth gape width in small fish of 15 mm standard length, decreasing to 50% of mouth gape width in fish of 40 mm standard length.